Low-Energy Electronic Structure of the High-T c 
Cuprates La2_ x Sr x Cu04 Studied by Angle-resolved 
Photoemission Spectroscopy 

T. Yoshida 1 , X. J. Zhou 2 , D. H. Lu 2 , Seiki Komiya 3 , Yoichi 
Ando 3 , H. Eisaki 4 , T. Kakeshita 5 , S. Uchida 6 , Z. Hussain 7 , 
Z.-X. Shen 2 and A. Fujimori 1 

department of Complexity Science and Engineering, University of Tokyo, Kashiwa, 
Chiba 277-8561, Japan 

2 Department of Applied Physics and Stanford Synchrotron Radiation Laboratory, 
Stanford University, Stanford, CA94305, U. S. A. 

3 Central Research Institute of Electric Power Industry, Komae, Tokyo 201-8511, 
Japan 

4 National Institute of Advanced Industrial Science and Technology, Tsukuba 
305-8568, Japan 

Superconductivity Research Laboratory, ISTEC, Shinonome, Koto-ku, Tokyo 
135-0062, Japan 

6 Department of Physics, University of Tokyo, Bunkyo-ku, Tokyo 113-0033, Japan 
7 Advanced Light Source, Lawrence Berkeley National Lab, Berkeley, CA 94720, U. S. 
A. 

E-mail: teppei@k .u-tokyo . ac . jp 

Abstract. We have performed a systematic angle- resolved photoemission spec- 
troscopy (ARPES) study of the high-T c cuprates La2- x Sr a; Cu04, ranging from the 
underdoped insulator to the superconductor to the overdoped metal. We have revealed 
a systematic doping evolution of the band dispersions and (underlying) Fermi surfaces, 
pseudogap and quasi-particle features under the influence of strong electron-electron 
interaction and electron-phonon interaction. The unusual transport and thermody- 
namic properties are explained by taking into account the pseudogap opening and the 
Fermi arc formation, due to which the carrier number decreases as the doped hole 
concentration decreases. 
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1. Introduction 

The mechanism that causes the high-temperature superconductivity in the copper oxide 
materials (cuprates) is still unknown. As the charge carriers (electrons or holes) are 
introduced into the parent antiferromagnetic insulator, the material evolves from an 
insulator to a superconductor, and eventually to a normal metal. Thus, the cuprate 
systems drastically change their behaviors depending on the carrier density in the two- 
dimensional CuC>2 planes, which are the stage of the high-T c superconductivity and 
related low-energy physics. While the overdoped (x > 0.2) cuprates behave like a 
relatively conventional Fermi liquid above T c , underdoped (x < 0.15) cuprates in the 
normal state show behaviors which are strongly deviated from those of a standard 
Fermi liquid. Therefore, it is necessary to know the detailed information about the 
doping evolution of the electronic structure of the CuC>2 plane to understand the cuprate 
systems. 

The recent remarkable development in the angle-resolved photoemission spec- 
troscopy (ARPES) technique has enabled us to investigate fine spectral features near 
the chemical potential, which may be used to interpret the thermodynamic and trans- 
port properties of solids from a microscopic point of view. Indeed, ARPES studies on 
cuprate superconductors, largely on the Bi 2 Sr 2 CaCu208 + j / (Bi2212) family, have been 
performed extensively pQ and revealed the normal-state electronic structure such as the 
dispersions of quasi-particle bands crossing the chemical potential, the presence of a 
pseudo-gap above T c and its d-wave like gap anisotropy in the underdoped region. On 
the other hand, bilayer splitting BiO-layer structural modulation, etc., complicate 
the interpretation of the ARPES spectra of Bi2212. 

In order to avoid the above complications and to study the intrinsic doping 
dependence of the electronic structure over a much wider hole concentration range, the 
La2- x Sr x Cu04 (LSCO) system is the most suitable one among the families of high-T c 
cuprate superconductors. LSCO has a simple crystal structure with single Cu0 2 layers, 
and no Cu-0 chains as in YBa2Cu307_j, (YBCO) nor complicated structural modulation 
of the block layers as in Bi2212. The hole concentration in the Cu02 plane can be 
controlled over a wide range and uniquely determined by the Sr concentration x (and the 
small oxygen non-stoichiometry) . One can therefore investigate the doping dependence 
from the undoped insulator (x = 0) to the heavily overdoped metal (x = 0.3) in a single 
system. Indeed, the most extensive set of data on the doping dependence of the transport 
and thermodynamic properties are available for the LSCO system [HI IH |KJ E] 
and can be in principle compared with ARPES data. 

ARPES studies of LSCO had been behind those of Bi2212 for some time due to 
the experimental difficulties in obtaining atomically flat surfaces by cleaving and to 
the relatively poor surface chemical stability. The first successful ARPES study of 
LSCO was performed by Ino et al. jTHj. They showed an overall feature of band 
dispersions and Fermi surface, including a crossover from a hole-like to an electron- 
like Fermi surface with hole doping [TQl [TTJ [TJ| . Peculiar features such as an unusual 
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suppression of the zone diagonal direction, which had not been observed in Bi2212, 
were also reported in [TU] and subsequently by Zhou et al. |13 3 - In addition, so-called 
"two-component" behavior, namely, the coexisting remnant lower Hubbard band and 
quasi-particle band with the chemical potential pinned within the charge-transfer gap, 
have been observed in the underdoped region ^Tj. Those unconventional behaviors were 
discussed in terms of microscopic charge inhomogeneity or the spin-charge fluctuations 
in a stripe form ^3JE3, which has been suggested from incommensurate peaks in the 
inelastic neutron scattering studies ^Hj- Alternatively, the apparent "two-component" 
behaviors in the ARPES line shapes have been interpreted by the strong electron-phonon 

coupling [muni eh]. 

Because of the weakness of the nodal spectral weight, concerns about the sample 
surface quality had persisted. Also, there were questions about the effects of 
photoemission matrix elements [20], which makes it difficult to extract quantitative 
information about the stripe effects from comparison between experiment and theoretical 
studies which predicted the suppression of spectral weight in the nodal direction 
[2111221123121]. However, recent ARPES studies on LSCO have revealed that under a 
suitable experimental geometry, a quasi-particle (QP) band with a clear kink feature 
crosses the Fermi level (Ep) in the nodal (0,0)-(7r,7r) direction [23] even down to the 
lightly doped region It has also been established that the Fermi surface is basically 
large [21] and that the (7r,0) region remains (pseudo) gapped, leading to the electronic 
state in which only part of the Fermi surface survives as an "arc" around the nodal 
direction. Such a large Fermi surface picture with the pseudogap has been theoretically 
studied both from RVB [2H] and Fermi liquid approach [21]. In this papar, we present a 
large set of ARPES spectra of LSCO showing the QP structure or the pseudogap feature 
on the entire Fermi surface, and compare the ARPES spectra with the transport and 
thermodynamic properties in the entire doping range. Particularly, we demonstrate 
that, in the underdoped region, both the ARPES spectra and the physical properties 
are highly affected by the opening of the pseudogap and the concomitant formation of 
the Fermi arc. 

Figure ^ summarizes the pseudogap behaviors seen in the thermodynamic and 
transport properties of LSCO. As shown in the figure, the electronic specific coefficient 7 
[H] and the Pauli paramagnetic suseptibility x decrease with decreasing x, in contrast 
to the mass enhancement toward the metal-insulator transition usually considered in the 
Fermi-liquid picture and experimentally seen in the Fermi-liquid system Lai_ x Sr x Ti03 
[3*U] . Furthermore, although the Fermi surface seen by ARPES is large, the carrier 
density n deduced from the Hall coefficient Rh [HI] and the superffuid density n s from 
magnetic-field penetration depth [32] measurements decrease with decreasing x. That 
is, the carrier number n appears to be given by the hole concentration x rather than the 
band filling 1 — x in the underdoped region |3*T] I33j . These unconventional behavior may 
be attributed to the pseudogap opening through which a portion of the Fermi surface 
becomes gapped above T c . The pseudogap opening provides a natural scenario to the 
long-standing issue of how the electronic structure evolves with hole doping from the 
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Figure 1. (Color online) Doping dependence of the electronic specific coefficient 7 UJ, 
Pauli paramagnetic suseptivility \ EH 1 t ne carrier density deduced from Hall coefficient 
Rh an d the superfluid density n s /m* [35] in LSCO. 

Mott insulator to the superconductor in two-dimensional systems. These features shall 
be described by the recent high-quality ARPES data of LSCO as described below. 

2. Experiment 

The ARPES measurements were carried out at beamline 10.0.01 of Advanced Light 
Source (ALS) and beamline 5-4 of Stanford Synchrotron Radiation Laboratory (SSRL), 
using incident photons with energies of 55.5 eV and 22.4 eV, respectively. SCIENTA 
SES-2002 (ALS) and SES-200 (SSRL) spectrometers were used in the angle mode. The 
total energy resolutions were about 15 meV (SSRL) to 20 meV (ALS). High-quality 
single crystals of LSCO were grown by the traveling-solvent floating-zone method. The 
samples were cleaved in situ and measurements were performed at about 20 K (ALS) 
or 10 K (SSRL). In the measurements at ALS, the electric field vector E of the incident 
photons lies in the Cu0 2 plane, rotated by 45 degrees from the Cu-0 bond direction, 
so that its direction is parallel to the Fermi surface segment around the nodal region. 
This measurement geometry enhances dipole matrix elements in this k region |2"H] . 
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Figure 2. Photoemission spectra of La2CuC>4 along the (0,0)-(7r,7r) nodal direction in 
the first Brillouin zone |18|. The corresponding momentum runs from (0.32,0.32)77 to 
(0.76,0.76)77, as indicated by the arrow, (a) Raw data, (b) To highlight the momentum 
dependence, a "background" given by a spectrum near (77, 7r) has been subtracted. In 
both (a) and (b), the bold curves correspond to a momentum closest to (77/2,77/2). 



3. Doping evolution of the electronic structure 

3.1. Doping Evolution of Quasi-particles 

One of the central issues in the field of high T c cuprates is how an antiferromagnetic 
insulator evolves with hole doping. Figure El reveals some important features such as 
the doping- dependent dispersions along the (remnant) Fermi surface and spectral weight 
transfer from the high-energy Hubbard band to the low-energy QP band. In order to 
further look into this problem, first, we show the energy distribution curves (EDC's) of 
the lowest energy feature along the nodal direction of undoped La 2 Cu0 4 . One can see 
one broad dispersive feature approximately folded back with respect to the (7r/2,7r/2) 
point [IS], similar to the results for Ca 2 Cu0 2 Cl 2 [IZj. As in the case of Sr 2 Cu0 2 Cl 2 , the 
dispersion of this feature agrees with that of t — t' — t" — J model well [SU EHH EH] • The 
broad band at ~ 0.5 eV has been interpreted as a boson side band showing a Franck- 
Condon broadening. A shell-model calculation showed that the width of the phonon 
sideband is in good agreement with the ARPES data of undoped La 2 Cu0 4 [IS]. This 
indicates that the electron-phonon coupling is strong enough, which seems a universal 
feature in the undoped cuprates. 

Next, let us look at the doping evolution of the spectral intensity of the QP in 
LSCO. FigureEfa) and (b) shows the evolution of spectra at k ~ (vr/2, 7r/2) and (vr, 0), 
respectively, with hole doping [2B]. At ~ (vr/2, vr/2), a finite spectral weight exists at E F 
except for x = 0, and increases with x without an abrupt change across the "insulator" - 
superconductor transition boundary of x ~ 0.06. Note that, even in the lightly-doped 
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Figure 3. (Color online) ARPES spectra at k = kp in the nodal direction in the 
second Brillouin zone (BZ) (a) and those at (tt, 0) (b) for various doping levels |26j. 
(c) Doping dependence of the nodal QP spectral weight, -Znqp- 



non-superconducting x=0.03, a QP peak is depressed and crosses the chemical potential 
in the nodal direction and forms a Fermi "arc" . This would explain the metallic behavior 
at high temperatures of the lightly-doped materials |37|. To see the doping evolution of 
QP, the spectral weight is plotted in Fig. |^c) as a function of hole doping x. Here, the 
nodal QP weight Z N q P is defined as the peak intensity of the EDC at k F in Fig.E^a). In 
the underdoped regime, -Znqp increase rapidly with hole doping, exibiting an apparent 
similarity to the doping evolution of the Drude weight in the optical conductivity [33J. 
The observed finite spectral weight at x=0.03 is consistent with the optical study of 
lightly-doped LSCO, which has shown that the Drude weight is finite already in x = 

0.03 t HH]. 

3.2. Fermi Surface and Two-dimensional Band Structure 

In order to deduce the Fermi surface or "underlying" Fermi surface, we have performed 
spectral-weight mapping at the chemical potential in the momentum space for various 
doping levels as shown in Fig. [T] [2?] . The spectral- weight distribution of the overdoped 
£=0.3 as well as 2=0.22 samples clearly shows an electron-like Fermi surface |25j . 
Spectral weight near (ir,0), a remnant (n,0) feature of the flat band becomes 
weaker in going from x=0.22 to £=0.30, indicating that, the more heavily the samples 
are overdoped, the more the electronic structure becomes conventional Fermi-liquid- 
like. The Fermi surface topology changes from electron-like in x=0.22 to hole-like in 
£=0.15, which may cause an electronic topological transition (ETT) |3E]. Below x=0.15, 
the nodal direction in the second Brillouin zone (BZ) become more intense than that 
around (7r,0), in contrast to x=0.22 and 0.3, due to the superconducting or pseudogap 
opening with d-wave symmetry in the underdoped and optimally doped regions. In the 
lightly-doped x=0.03 [panel (a)], the (ir,0) spectral weight almost disappears, and only 
the spectral weight around the zone diagonal direction remains appreciable around the 
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chemical potential. 

The EDO's along the (0, 0) — (jr, 0) — (jr, ti) lines for various doping levels are shown 
in Fig. Eland their intensities along the same lines are plotted in the energy- momentum 
E — k space in Fig. |3 One can see that the flat band feature appears around (ir, 0) 
in the optimally and underdoped region x < 0.15, while the QP dispersion crosses the 
chemical potential along the (0,0)-(7r, ti) line. Also, the intensity plots clearly show how 
the intensity at the chemical potential around (vr,0) becomes weak for x < 0.15 as the 
pseudo-gap and/or the superconducting gap opens around (n,0). One can see from 
Fig. El that the line width around (vr,0) becomes broad with decreasing x, indicating 
increasing antiferromagnetic correlation with decreasing x. Also, the "two-component" 
behavior is clearly seen in the EDO's of the lightly-doped samples with £=0.03 and 0.05. 

The "two-component" behavior in LSCO ^lj EHj was interpreted in terms of 
microscopic charge inhomogeneity or stripe. According to a recent theoretical study, 
a mixed phase of antiferromagnetic and superconducting states has been proposed to 
explain the characteristic "two-component" behavior in LSCO j3Hl- At the same time, a 
new way to interpret the "two- component" spectral feature in terms of polaronic effects 
has also emerged fTf \ I19 [ IT%]. In the new picture, the high energy ~ 0.5 eV have been 
interpreted as polaronic side bands. These two interpretations may be related to each 
other because it has been predicted that phase separation between the insulating and 
metallic phases occurs under a certain regime of electron-phonon coupling strengths 

SB- 

The Fermi surfaces obtained from the two-dimensional tight-binding fits are 
superimposed on the spectral weight mapping in Fig. The experimental hp's have 
been fitted to the two-dimensional (2D) single-band tight-binding (TB) model 

Ek = — 2t[cos(k x a) + cosf^a)] — 4t' cos(k x a) cos(k y a) 

- 2t"[cos(2k x a) + cos(2k y a)} + e , (1) 

as shown by blue curves. Here, t, t' and t" are the first, second and third nearest 
neighbor transfer integrals between Cu sites. We have assumed constant t = 0.25 eV and 
relationship —if' /if =1/2 for all the doping levels, and regarded t' and Eq as adjustable 
parameters. The Fermi surface could be well fitted by the TB model, although some 
misfit can be seen for the kink structure, which arises from coupling to phonons as 
discussed in Sec. 4 [121 HH] and the extremely flat band dispersion around (7r,0) in the 
underdoped region as shown in FigH3 In Fig. E]T), the deduced Fermi surfaces for 
various x's are replotted, clearly showing the change of the Fermi surface shape and 
topology with x. The TB dispersions with the same parameters are superimposed on 
the color mapping in Fig. |H1 Overall features of the Fermi surfaces and band dispersions 
are in good accordance with the results of the local-density approximation (LDA) band- 
structure calculations |44j . 

The obtained TB parameters —t'/t shown in Fig. 0] exhibit a clear doping 
dependence. The increase of —t'/t with decreasing x can be qualitatively explained 
by the increase of the Cu-apical oxygen distance with decreasing x [15] • In FigJZ^a), 
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Figure 4. (Color) Spectral weight mapping in fc-space at the chemical potential in 
LSCO. Data were taken with hv=hh.h eV. Red dots indicate kp positions determined 
by the peaks of momentum distribution curves (MDC's) at the chemical potential. 
The blue curves show the Fermi surface interpolated using the tight-binding model 
|27|. The evolution of the Fermi surface with doping is more clearly seen in panel (f). 

the —t'/t values are plotted as a function of the Cu-apical oxygen distance J23] and 
compared with those of the LDA calculation by Pavarini et al. [46 . The obtained 
—t'/t agree rather well with the results of the LDA calculations —t'/t ~ 0.15 [46J. 
The LDA calculation has indicated that the larger Cu-apical oxygen causes the larger 
—t'/t, which is consistent with the trend in the obtained —t'/t values deduced here. The 
present result is also consistent with the cellular dynamical mean field theory (CDMFT) 
calculation ^J\, in which the momentum dependence of the self-energy has been taken 
into account by using a cluster instead of a single-impurity atom in DMFT. As shown 
in Fig. Ulb), the CDMFT calculation has indicated that increasing correlation effects 
in the underdoped region is reflected on the effective increase of the —t'/t value. It 
therefore seems that both the structural effects as predicted by the LDA calculation 
and the correlation effects as predicted by the CDMFT calculation increase —t'/t as x 
decreases and explain the ARPES results. 

As for Luttinger's sum rule, the underlying "Fermi surface" of LSCO determined 
from the low energy spectral weight of ARPES approximately fulfills it remarkably well 
even down to the lightly-doped region as shown in Fig. Et[a) [27 . This is in strong 
contrast with the results on Na-CCOC, which shows a clear deviation from Luttinger's 
sum rule in the underdoped region |4"8] . These differences between LSCO and Na-CCOC 
may be related to the difference in the t'/t value and/or to the difference in the electron- 
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Energy relative to E F (eV) 

Figure 5. Energy distribution curves (EDC's) along the (0,Q)-(7r, 0)-(ir, ir) lines for 
various doping levels in LSCO. Thick curves represent the spectra at (ir, 0). The spectra 
along the (0,0)-(7r,0) line are from the second BZ. 

lattice coupling strength. In relation to the area of the Fermi surface, the shape of the 
Fermi surface may be characterized by the position of the Fermi surface along the nodal 
direction. As shown in Fig. |Hl(b), the same position in Na-CCOC shows a strong doping 
dependence and is extrapolated to (7r/2,7r/2) in the x=0 limit, whereas in LSCO, the 
limiting value is off the (tt/2,tt/2) point ^T7\. This may be related to the observation that 
the chemical potential approaches the top of the lower Hubbard band in Na-CCOC as 
the doping decreases whereas the chemical potential stays away from the lower Hubbard 
band in LSCO HHUnHHD!. 

3.3. Large Pseudogap and Small Pseudogap 

There are two kinds of "pseudogaps" in the high-T c cuprates called a "small pseudogap" 
and a "large pseudogap". First, we discuss the relationship between the high-energy 
pseudogap behaviors, which have been observed as a characteristic temperature in the 
magnetic susceptibility x(T) [H] and the Hall coefficient Rh(T) and a characteristic 
energy in the photoemission spectra. As shown in Fig. Efb), the doping dependence of 
the flat band position -E^o) agrees with that of T x below which x(T) decreases from the 
maximum value, if the relationship for the ci-wave mean-field relation Eu^\ = 4.3fceT/2 
is assumed. Figure EJa) shows angle-integrated spectra obtained by integrating ARPES 
spectra in the second BZ, where the -E^o) is indicated by vertical bars. Since the flat 
band gives a singularity in the density of states (DOS), the corresponding maximum 
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Figure 6. (Color) Intensity plot in E-k space along the symmetry lines (0,Q)-(tt,0)- 
(ir, 7r)-(0,0) in LSCO. Data were taken with hv=hh.b eV. Blue curves show the tight- 
binding band dispersion with the parameters shown in Fig. 



is seen at -E^o) in the angle-integrated spectra j^U This energy scale may be 

related to the strong suppression of the spectral weight over an extended energy range, 
particularly prominent near (tt,0). These results are consistent with the theoretical 
prediction of ETT that E( n ^ is a crucial factor to determine the quantum critical point 

Next, we discuss the relationship between the small pseudo gap behaviors of the 
thermodynamic and transport properties and the low-energy spectral feature of ARPES 
spectra. Figure ITOT a) shows ARPES spectra on the Fermi surface in the antinodal 
direction. In order to clearly show the gap opening, the spectra in Fig. ITUT a) have 
been symmetrized with respect to Ep as shown in Fig. fTDT bh For x= 0.3 and 0.22, 
the symmetrized spectra show a peak at the chemical potential in contrast to those 
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Figure 7. (Color online) Comparison of the —t'/t values deduced from the tight- 
binding fit to the ARPES data with theoretical predictions, (a) Comparison between 
—t'/t calculated by LDA 46 (filled diamond) and those estimated from ARPES (circle), 
plotted as a function of the Cu-apical oxygen distance The largest and smallest 
—t'/t in the ARPES data corresponds to those of a;=0.03 and £=0.30 samples, 
respectively, (b) Comparison of t'/t estimated from ARPES Fermi surface and those 
calculated using CDMFT 02] ■ The theoretical values have been normalized to the bare 
(i.e., input) t'/t, while the experimental values have been normalized to the values at 
£=0.3. n — I — x is the electron number. 
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Figure 8. (Color online) Fermi surface area and shape of LSCO as functions of 
doping |27]. Data for Na-CCOC @H] are also plotted, (a) Doping dependence of the 
hole number xps deduced from the Fermi surface area. Luttinger's sum rule xfs = x 
is shown for comparison, (b) Fermi surface position in the nodal direction. 
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Figure 9. Angle-integrated photoemission spectra and the large pseudogap. (a) 
Spectrum for each doping. Vertical bars represent the binding energy of the flat band, 
(b) Doping dependence of the binding energy of the flat band is compared with the 
characteristic temperature T x where the magnetic susceptibility x(T) takes a maximum 

m 



for x=0.15, corresponding to the normal Fermi-liquid behavior of the transport results 
[3T] . Figure IToT c) shows the doping dependence of the relative energy of the leading 
edge midpoint (LEM). Although the LEM in x=0.22 is above the Fermi level in spite of 
T c > 10K, it is by ~ (2 — 3) meV below that of the non-superconducting sample £=0.3, 
implying the existence of a small gap. Note that, when the QP is at the Fermi level, 
the LEM would be above the Fermi level due to the energy broadening under the finite 
energy resolution. As seen in Fig. HUT c). the observed gap well corresponds to the small 
pseudogap observed by the scanning tunneling measurements. 

4. Effects of Electron-Phonon Coupling 

It has been recognized from the very beginning of the high-T c research that many- 
body effects are key to understanding cuprate physics. Due to its proximity to the 
antiferromagnetic Mott insulating state, electron-electron interactions are extensively 
discussed in the literaturepjj EH]. In this Section, we shall briefly describe the 
present status of our understanding of electrons interacting with phonons. A more 
comprehensive review is presented elsewhere [56;. 

4-1. Electron-Phonon coupling along the Nodal Direction 

The d-wave superconducting gap is zero along the nodal direction. As shown in 
Fig lTlT a). the energy- momentum dispersion curves from MDC method exhibit an abrupt 
slope change ( "kink" ) near 70 meV. The kink is accompanied by an accelerated drop in 
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Figure 10. ARPES spectra on the Fermi surface in the antinodal region, (a) EDC 
for each doping, (b) Symmetrized spectra of (a) . (c) Relative shift of the leading-edge 
midpoint (LEM) in the ARPES spectra in (a). The LEM shift is compared with the T c 
(filled square) and the gap measured by scanning tunnel spectroscopy (filled triangles) 
|53j . assuming the relationship E = A.SkgT between energy and temperature 54 . 

the MDC width at a similar energy scale (Fig. fTTT b)). The existence of the kink has been 
well established as ubiquitous in hole-doped cuprate materials [5T | W2 \ 15% } IHO j IfiT] : 

1. It is present in various hole-doped cuprate materials, including Bi2212, 
Bi 2 Sr 2 Cu06 (Bi2201), LSCO and others. The energy scale (in the range of 50-70 meV) 
at which the kink occurs is similar for various systems. 

2. It is present both below T c and above T c . 

3. It is present over an entire doping range (Fig. ITTT a)). The kink effect is stronger 
in the underdoped region and gets weaker with increasing doping. 

While there is a consensus on the data, the exact meaning of the data is still under 
discussion. The first issue concerns whether the kink in the normal state is related 
to an energy scale. Valla et al. [62J argued that the system is quantum critical and 
thus has no energy scale, even though a band renormalization is present in the data. 
Since their data do not show a sudden change in the scattering rate at the corresponding 
energy, they attributed the kink in Bi2212 above T c to the marginal Fermi liquid (MFL) 
behavior without an energy scale |59j. Others believe the existence of energy scale in the 
normal and superconducting states has a common origin, i.e., coupling of quasiparticles 
with low-energy collective excitations (bosons) H21 EH] • The sharp kink structure in 
dispersion and concomitant existence of a drop in the scattering rate which is becoming 
increasingly clear with the improvement of signal to noise in the data, as exemplified in 
underdoped LSCO (x=0.063) in the normal state (Fig. ITTTb))[13|. are apparently hard 
to reconcile with the MFL behavior. 

A further issue concerns the origin of the bosons involved in the coupling, with a 
magnetic resonance mode [5E1 EHj and optical phonons [12] being possible candidates 
considered. The phonon interpretation is based on the fact that the sudden band 
renormalization (or "kink") effect is seen for different cuprate materials, at different 
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Figure 11. (Color online) Doping dependence of the nodal QP in LSCO 43 . (a) 
Dispersion of LSCO with various doping levels (x=0.03 to 0.30) measured at 20 K 
along the (0,0)- (71771") nodal direction, (b). Scattering rate as measured by MDC width 
(full-width-at-half-maximum, FWHM) of the LSCO (x = 0.063) sample measured at 
20 K. 

temperatures, and at a similar energy scale over the entire doping range j32]- For 
the nodal kink, the phonon considered in the early work was the half-breathing mode, 
which shows an anomaly in neutron experiments jOHUM!- Unlike the phonons, which are 
similar in all cuprates, the magnetic resonance (at correct energy) is observed only in 
certain materials and only below T c . The absence of the magnetic mode in LSCO and 
the appearance of the magnetic mode only below T c in some cuprate materials are not 
consistent with its being the cause of the universal presence of the kink effect. Whether 
the magnetic resonance can cause any additional effect is still an active research topic 

While the nodal data clearly reveal the presence of coupling to collective modes 
with well-defined energy scale, there are a couple of peculiar behaviors associated 
with the doping evolution of the nodal dispersion. As seen from Fig. ITlT a). the 
low-energy dispersion and velocity from the MDC method is insensitive to doping 
over the entire doping range, while the high energy velocity increases with decreasing 
doping. [13] • In terms of conventional electron-phonon coupling, if one considers that the 
"bare band" does not change with doping but the electron-phonon coupling strength 
increases with decreasing doping, as it is probably the case for LSCO, one would 
expect that the low energy dispersion and velocity show strong doping dependence, 
while the high-energy ones converge. This is opposite to the experimental observations. 
Moreover, considering the electron-electron interaction gets stronger with decreasing 
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doping, according to conventional wisdom, this would result in a larger effective mass and 
smaller velocity, which is again opposite to experimental observation (Fig. ITlT a)). These 
anomalies indicate a potential deviation from the standard Migdal-Eliashberg theory 
and the possibility of a complex interplay between electron-electron and electron-phonon 
interactions. This phenomenon is a hint of polaronic effect where the traditional analysis 
fails. Such a polaron effect gets stronger in deeply underdoped system even along the 
nodal direction. It is noted that, while the low energy velocity is rather "universal" 
via MDC analysis, the low energy dispersions derived from EDC and MDC methods 
are rather different JH]- The difference is most pronounced at low doping where the 
coupling is strong. In this case, the EDC-derived velocity exhibits considerable doping 
dependence. While for weakly interacting systems MDC and EDC are expected to 
yield similar dispersion relations from linearly dispersing bands, strong coupling features 
which cause substantial renomaliations may cloud EDC and MDC analysis. The origin 
of the difference between MDC and EDC-derived dispersions remains a matter to be 
further explored. 

4-2. Multiple Modes in the Electron Self-Energy 

In conventional superconductors, the successful extraction of the phonon spectral 
function, or the Eliashberg function, a 2 F(u), from electron tunneling data played 
a decisive role in cementing the consensus on the phonon-mediated mechanism of 
superconductivity [S7| . For high temperature superconductors, the extraction of the 
bosonic spectral function can provide fingerprints for more definitive identification of 
the nature of bosons involved in the coupling. 

In principle, the ability to directly measure the dispersion, and therefore, the 
electron self-energy, would make ARPES the most direct way of extracting the bosonic 
spectral function. This is because, in metals, the real part of the electron self-energy 
ReS is related to the Eliashberg function a 2 F(fl; e, k) by: 

ReS(k, e; T) = ^ dQa 2 F(Q; e, k)K ( — , — J , (2) 

where 

/oo 2z 
dx— -f(x + y), (3) 
-00 x 2 — z l 

with f[x) being the Fermi distribution function. Such a relation can be extended 
to any electron-boson coupling system and the function a 2 F(u) then describes the 
underlying bosonic spectral function. Unfortunately, given that the experimental data 
inevitably have noise, the traditional least-square method to invert an integral problem 
is mathematically unstable. 

Very recently, Shi et al. have made an important advance in extracting the 
Eliashberg function from ARPES data by employing the maximum entropy method 
(MEM) and successfully applied the method to Be surface states |68j . The MEM 
approach jHEj is advantageous over the least squares method in that: (i) It treats the 
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bosonic spectral function to be extracted as a probability function and tries to obtain 
the most probable one. (ii) More importantly, it is a natural way to incorporate the 
priori knowledge as a constraint into the fitting process. In practice, to achieve an 
unbiased interpretation of data, only a few basic physical constraints to the bosonic 
spectral function are imposed: (A) It is positive. (B) It vanishes at the limit lj^O. (C) 
It vanishes above the maximum energy of the self-energy features. As shown in the case 
of Be surface state, this method is robust in extracting the Eliasberg function [HS]. 

Initial efforts have been made to extend this approach to underdoped LSCO and 
evidence for electron coupling to several phonon modes has been revealed [HH|. As seen 
from Fig. from both the electron self-energy (Fig. IT27 a)). and the derivative of their 
fitted curves ((Fig. IT2T a)). one can identify two dominant features near ~ 40 meV 
and ~60 meV. In addition, two addition modes may also be present near ~25meV and 
~75 meV [OH]- The multiple features in Fig. IT27 b) show marked difference from the 
magnetic excitation spectra measured in LSCO which is mostly featureless and doping 
dependent [ZD]- In comparison, they show more resemblance to the phonon density-of- 
states (DOS), measured from neutron scattering on LSCO (FigJT^c)) [7T], in the sense 
of the number of modes and their positions. This similarity between the extracted fine 
structure and the measured phonon features favors phonons as the bosons coupling to 
the electrons. In this case, in addition to the half-breathing mode at 70~80 meV that 
we previously considered strongly coupled to electrons [12], the present results suggest 
that several lower energy optical phonons of oxygens are also actively involved. 

We note that, in order to be able to identify fine structure in the electron self-energy, 
it is imperative to have both high energy resolution and high statistics |72j . These 
requirements have made the experiment highly challenging because of the necessity 
to compromise between two conflicting requirements for the synchrotron light source: 
high energy resolution and high photon flux. Further improvements in photoemission 
experiments will likely enable a detailed understanding of the boson modes coupled to 
electrons, and provide critical information for the pairing mechanism. 



5. Relation to thermodynamic and transport properties 

The thermodynamics and transport properties of the high-T c cuprates are fundamentally 
affected by the strong electron correlation near the filling-control Mott transition, 
particularly by the opening of the pseudogap as implied by FigJT] The electronic specific 
heat coefficient 7 decreases with decreasing hole concentration in the underdoped region 
[5|. The electrical resistivity shows unconventional electrical resistivity such as the 
T-liner temperature dependence in the optimally doped region [21] and the metallic 
behavior with resistivity well exceeding the Ioffe-Regel limit in the lightly doped region 
|37j . In order to elucidate the origin of these unconventional phenomena, in this 
Section, we shall attempt to explain the doping dependence of the electronic specific 
heat coefficient 7, superfluid density, the electrical resistivity and the Hall coefficient on 
the basis of the ARPES results. 
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Figure 12. (Color online) Multiple modes coupling in the electron self-energy in 
LSCO 6.9;- (a) Effective real part of the electron self-energy for LSCO x=0.03 (al), 
0.063 (a2), 0.07 (a3)and ^0.06 (a4) samples, (b). Second-order derivative of the 
calculated ReE. The four shaded areas correspond to energies of 23-29, 40-46, 58-63 
and 75-85 meV where the fine features fall in. (c) Phonon density of states F(us) for 
LSCO x = (red) and x = 0.08 (blue) measured by neutron scattering |71j . 




Figure 13. (Color online) Symmetrized EDC's along the (underlying) Fermi surface 
for x=Q. 03(a), x=Q. 063(b) and x=Q. 15(c). Panel (d) shows the arc length ratio 
defined by the fraction of the Fermi surface length which shows a QP peak as in 
the symmetrized EDC's. Curves for n = x and nnaii = 1/gRh |31| are also shown. 
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Figure 14. (Color online) Density of states at the chemical potential in LSCO. (a) 
Fermi velocity vp along the Fermi surface, (b) Comparison the electronic specific heat 
coefficient 7 [5] and the QP DOS at the chemical potential calculated using vf from 
ARPES data. The QP DOS from ARPES with the arc length taken into account is 
also shown. 



5.1. Electronic specific heat 

The electronic specific heat coefficient 7 of a Fermi liquid is given by 7 = (n 2 k%/3)N(0)* , 
where N(0)* is the QP density of states (DOS) at the chemical potential. In a single- 
band system like single-layer cuprates, each momentum k contributes one QP, if the 
system is a Fermi liquid. Therefore, by using the Fermi velocity Vp of the QP band 
along the Fermi surface, one can deduce the two-dimensional QP DOS iV(0)* using 
the formula N(0)* = (1/2tc 2 ) J ds/hvp, where the integrations is made over the Fermi 
surface. 

In order to take into account the effect of the pseudogap in the present analysis, 
we first attempted to derive the length of the Fermi arc. From symmetrized EDO's on 
the (underlying) Fermi surface as shown in Fig. EH ( a )-( c )? we have defined the arc 
length by the portion of the (underlying) Fermi surface showing a single QP peak at the 
chemical potential. As shown in Fig. iTHT d). one can see that the arc length evolves with 
hole doping systematically and eventually it becomes a fully extended (100%) closed 
Fermi surface in overdoped x=0.22. The carrier number defined by nn a ii = 1/eRn are 
compared with the arc length in Fig. iTHT d). The difference between r^aii and the arc 
length may come from the fact that nnaii is affected by the curvature of the Fermi 
surface, too [ZH]- Note, however, that due to the finite energy resolution, closely spaced 
two QP peaks may not be resolved in the symmetrized EDO's and that the arc length 
may be overestimated. 

Figure fTlT a) shows the Fermi velocity vp for various doping levels determined from 
the observed QP dispersions as a function of Fermi-surface angle a. Here, we have used 
MDC for determining vp in a similar manner as in Fig. ^2 while an alternative way is to 
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use EDC's [71]. One can see that the vp remains relatively unchanged with hole doping 
not only in the node direction (a = 45°; see Fig. ITTj) but also on the entire Fermi surface. 
In Fig. IT4T b). the electronic specific heat coefficient 7 calculated using the Vp values 
thus deduced is compared with the 7 of LSCO samples where the T c was suppressed 
by Zn-doping [*"]. The iV(0)* calculated from ARPES has a maximum around x—0.15 
because the flat band near (tt, 0) becomes closest to the Fermi level. If the system is 
a Fermi liquid, N(0)* from ARPES and 7 should agree with each other. However, 7 
decreases much faster than N(0)* from ARPES with decreasing x, indicating that the 
QP density at the chemical potential is indeed depleted in the underdoped materials and 
demonstrating that the opening of the pseudogap around (tt, 0) affects the low-energy 
QP excitations. If the shrinking arc length is taken into account, one can qualitatively 
reproduce the decrease of 7 with decreasing x. As for the absolute values, however, the 
experimental values are larger than those deduced from ARPES by a factor of ~2. The 
origin of this discrepancy is not known at present. 

5.2. Superfluid density 

Comparison between the ARPES spectra and the superfluid density would also 
give important information about the pseudogap opening and its effects on the 
superconductivity "7J"]. The in-plane London penetration depth \ a b and the superfluid 
density p s = n s /m* are related through l/\ 2 ab = Hoe 2 n s /m* , where /x is the vacuum 
permeability, and n s is the density of superconducting electrons. For the density of 
superconducting electrons, we assumed that n s per Cu atom is given by (1 — a;)x(arc 
length), where the doping dependence of the arc length is given in Fig fTBT d). Since 
the nodal direction rather than the antinodal direction is expected to be the main 
contribution to the in-plane superfluid density due to the lighter effective mass, m* was 
approximated by m* = hkp/vF, where Vp is the Fermi velocity in the nodal direction 
[4*3] (Sec. 4.1). In Fig. [""""] the doping dependence of 1/A^, [7"""| is compared with that 
deduced from the ARPES data. The two curves show qualitatively good agreement 
between them. Considering that m* is nearly independent of doping [43 , the doping 
dependence of the arc length is the principal origin of the doping dependence of n s . 
Quantitatively, however, the 1/A^ 6 values show a difference of factor ~5, implying 
that the effective density of superconducting electrons is reduced compared with the 
nominal carrier number even in the overdoped region. Also, note that the superfluid 
density becomes finite only for x > 0.05, while that estimated from ARPES shows finite 
value even in £=0.03. Another viewpoint to explain the reduced superfluid density 
is to attribute it to the reduced QP spectral weight Z rather than the enhanced m* 
"23 EI], which also explain the doping dependence of 1/A„ 6 [see Fig|3[c)]. The overall 
quantitative difference may come from the effect of Z. 
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Figure 15. (Color online) In-plane superfluid density deduced from the magnetic 
penetration depth X a b [ZSj compared with that estimated from the ARPES data. 



5.3. Electrical resistivity 

Because of the pseudo-gap opening around (vr,0), the electronic states around the node 
mainly contribute to the in-plane transport properties in the underdoped region. Here, 
we have attempted to estimate the in-plane electrical resistivity from the observed 
MDC width A A; = 1/4 [Fig. IToT a)]. where Ik is the mean-free-path, and the length 
of the "Fermi arc" obtained in Sec. 5.1 as follows. In the simple Drude formula, the 
resistivity is given by p = m* /ne 2 r = hkpAk/ne 2 . It is often assumed that n = x based 
on the Hall coefficient results [3T]. In the Boltzmann equation in two-dimension, the 
diagonal component of the conductivity is given by a xx = (e 2 /hn) j Ik cos 2 O^ds, where 
the integration is made over the Fermi surface, and 9k is the angle between = vp-r^ 
and the electric field taken parallel to the x-axis [Z3]- Using ARPES data of high-T c 
cuprates, the electrical resistivity at high temperatures, where the pseudogap disappears, 
has been calculated using this method |78j . In addition to this treatment, we have taken 
into account the effect of "Fermi arc" by reducing the integration length on the Fermi 
surface. 

The MDC width at the chemical potential is plotted in Fig. UHl a) as a function of 
Fermi-surface angle a. One can see that the width around the node (a=45°) depends 
on doping only weakly, which is consistent with the weak doping dependence of the 
charge mobility p = er/m* = e/hkpAk j^Zj. Using the data in Fig. IToT a). the doping 
dependence of the in-plane electrical resistivity p was estimated by applying the Drude 
formula with n = x and the nodal MDC width, and the Boltzmann formula with and 
without the Fermi arc. Then, they are compared with the transport results [37] in Fig. 
Mh)- 

In all the p estimated from ARPES, the discrepancy from the transport data is 
pronounced in the over doped region. The discrepancy stems from the wide MDC 
widths compared to that expected from the transport results in the overdoped region, 
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Figure 16. (Color online) Electrical resistivity of LSCO estimated from ARPES. 
(a) MDC width for various doping levels plotted as a function of Fermi-surface 
angle a. (b) Comparison of the resistivity estimated using the MDC width and the 
electrical resistivity at 20 K 37 . For the transport data for those compositions which 
show localization behavior at low temperatures, the minimum resistivity at a higher 
temperature is adopted. 

which may be due to the effects of vertex correction, the possible dominance of 
forward scattering j75] in the overdoped region and/or the finite instrumental resolution 
neglected in the present analysis. On the other hand, in the underdoped region, the p 
estimated from ARPES including the effects of "Fermi arc" or assuming n = x is in 
good agreement with the transport results, while that evaluated using the entire Fermi 
surface is by more than one order of magnitude smaller than the transport data. This 
suggests that the pseudogap is the origin of the small carrier number and hence the high 
electrical resistivity of the underdoped samples. As a whole, the p values estimated from 
ARPES with the Fermi arc well explains the doping dependence of the transport data 
of the underdoped LSCO. 

5.4- Incommensurability in neutron scattering data 

Incommensurate peaks seen in the neutron scattering experiments of LSCO ^H] have 
been discussed as a signature of the spin-charge fluctuations in a stripe form [TH UK]. 
There have been several attempts to explain the incommensurability from the shape of 
the Fermi surface l%H l%2"] . Recently, the QP interference has been observed in the 
scanning tunneling microscopy (STM) study of Bi2212 [SHI Ell- in another STM study, 
the observed features were discussed in connection with the stripe feature |85j . The QP 
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interference has been explained by the autocorrelations of the ARPES intensity map in 
momentum space (HHl EI] ■ 

Therefore, we applied the autocorrelation analysis to the present ARPES results 
of LSCO and found that neutron scattering peak intensity Inr\;o(q, oj = 0) can be 
approximately expressed by the autocorrelation formula Imx (q, 0) ~ SkA(k, 0)A(k + 
q, 0) as shown in Fig. El Although this formula may not be sufficient for correlated 
electron systems, autocorrelation analysis may give some insight into the problem of 
the incommensurability. Figure El (b) shows that the obtained incommensurability 
increases with hole doping, although there are quantitative deviation from the neutron 
results. Also, it should be noted that neutron scattering for x=0.03 shows diagonal peaks 
whereas our incommensurate peaks remained vertical. In the present autocorrelation 
model, such an increase of the incommensurability can be attributed to the decreasing 
Fermi momentum /jp measured from the T point with hole doping [Fig. Efb)], since the 
incommensurability is given by the momentum vector connecting nearly parallel part of 
the Fermi surface around the node, as shown in the Fig. H7T c). The antinodal region does 
not contribute to the autocorrelation function because spectral weight in that region is 
very small due to the pseudogap opening. It is interesting to note that the 4a x 4a 
ordering observed in Na-CCOC [EE] m &y be related to the nesting vector q = ir/2a 
connecting the Fermi surface near (tt,0) |4~5| . The calculated peak width are broadened 
in the overdoped x=0.22, which is also consistent with the neutron results. The present 
analysis suggests that the shape of the Fermi surface of LSCO with small —t'/t may 
have a nesting instability of the observed incommensurability. On the quantitative level, 
however, the obtained incommensurability is a little different from the stripe picture, 
and therefore further analysis with more sophisticated models will be needed to address 
the stripe issue. 

6. Conclusion 

We have performed a systematic ARPES study of the high-T c cuprates LSCO from the 
undoped insulator to the overdoped metal. The doping dependence of the (underlying) 
Fermi surfaces, the QP band dispersion and the pseudogap behaviors were studied in 
detail. We observed a doping dependence of the TB parameter —t'/t, which is consistent 
with the LDA calculation, which takes into account the effect of apical oxygen, as well as 
with the CDMFT calculation, which takes into account the increased electron correlation 
with decreasing hole concentration. Also, the Fermi surface deduced using MDC's was 
found to approximately fulfill Luttinger's sum rule even down to the lightly-doped region 

As for the QP spectral feature, in the underdoped region, apparent two-component 
behaviors occurs with the chemical potential pinning [!§] • The Fermi velocity was nearly 
doping independent [13] and the spectral weight of the nodal QP increased with hole 
doping. The length of the Fermi arc increased with hole doping. The relationship 
between the large pseudo gap and the singularity in the DOS caused by the flat band 
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Figure 17. (Color online) Autocorrelation mapping of the ARPES data for LSCO. (a) 
Intensity plots of T,kA(k,0)A(k + q,0) for various doping, (b) Doping dependence of 
the incommensurability S obtained from the autocorrelation mapping, (c) Momentum 
vectors corresponding to the incommensurate peaks for a;=0.15. 



feature is discussed. The small pseudo gap observed around (tt,0) was found to be 
consistent with other experimental data. 

Finally, by using the information from the ARPES results, we have successfully 
reproduced the qualitative behavior of the thermodynamic and transport properties 
of LSCO. In the analysis, we demonstrated that the pseudogap and the Fermi arc 
feature observed by ARPES can explain the unconventional doping dependence of the 
thermodynamic and transport properties of LSCO. 

Although the analysis using the ARPES data can qualitatively explain many 
physical properties of the underdoped LSCO, quantitative discrepancies remain to be 
explained. Obviously the definition of the Fermi arc and its length is rather loose and 
only approximate at the present stage, and will have to be further examined on a firm 
ground. Also, the origin of the difference between of LSCO and other cuprate families 
such as Na-CCOC remains to be understood. The effect of different electron-phonon 
coupling strengths and the effect of apical oxygen leading to the different t'/t and t"/t 
values have to be investigated in future studies. 
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